Is the Mott transition relevant to f-electron metals ? 
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We study how a finite hybridization between a narrow correlated band and a wide conduction 
band affects the Mott transition. At zero temperature, the hybridization is found to be a relevant 
perturbation, so that the Mott transition is suppressed by Kondo screening. In contrast, a first- 
order transition remains at finite temperature, separating a local moment phase and a Kondo- 
screened phase. The first-order transition line terminates in two critical endpoints. Implications for 
experiments on f-electron materials such as the Cerium alloy Ceo.gLao.iTho.i are discussed. 
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The Mott transition, i.e how electrons evolve from lo- 
calized to itinerant as a function of an external parameter 
such as pressure, is a fundamental problem in condensed 
matter physics. It is a key phenomenon in d-electron ma- 
terials, such as transition-metal oxides, in which a set of 
bands with rf-character is well separated and close to the 
Fermi level. In this case, the phase with localized elec- 
trons is insulating, and that with itinerant electrons is 
metallic. Dynamical mean- field (DMFT) studies Q have 
deepened our understanding of this phenomenon, and led 
to many interesting experimental predictions which have 
been recently verified experimentally [2|. 

In several /-electron materials, a transition between a 
phase where /-electrons are more localized and another 
in which they are more itinerant is also observed (such as 
the isostructural 7- a phase transition of Cerium) 0, . 
In these materials however, there is a broad band with 
spd-character close to the Fermi level (in addition to the 
/-orbitals), and both phases are metallic. It was sug- 
gested early on by B. Johansson |]| that the concept of 
a Mott transition within the /-electron subspace may 
still be relevant in this context. A different view is the 
Kondo volume collapse (KVC) model [H[zl@], m which 
the transition is driven by the change in the hybridisa- 
tion between the two phases with different unit-cell vol- 
umes. In this picture, the broad band of conduction elec- 
trons plays a key role, while it is merely a spectator in 
the Mott picture. The interplay of these two mechan- 
ims in Cerium have recently received a great deal of at- 
tention, since DMFT provides a comprehensive frame- 
work in which these problems can be studied systemati- 

cany sesiueli. 

In this article, we study the localization-delocalization 
transition within a simple model, which nevertheless re- 
tains the key ingredients present in /-electron materials. 
The model interpolates between a Hubbard model for 
the /-orbital and the periodic Anderson model (PAM) in 
which this orbital is hybridized to a broad band. Our goal 
is to understand whether, in a purely electronic model, 
the Mott transition present when the /-band is isolated 
remains a robust feature in the presence of a finite hy- 



bridization to a broad band. Our key finding is that 
the answer to this question depends on temperature in a 
crucial manner. At zero-temperature, the Kondo effect 
always sets in and screens the local moment. As a result, 
in a purely electronic setting, the Mott transition is sup- 
pressed by an arbitrarily small hybridization. In contrast, 
a first- order transition remains at finite temperature. 

The similarities between the phase diagram of the 
PAM and that of the Hubbard model at finite tempera- 
ture have been pointed out in previous work |l3| . How- 
ever, the distinction between a first-order transition with 
coexisting electronic phases and a mere crossover was 
not addressed. More importantly, the zero-temperature 
case was investigated in the case where the hybridisa- 
tion vanishes at the Fermi level [Til fl5| : this is a non- 
generic case in which the transition survives down to 
T = 0. In the generic case of a finite hybridisation, 
the connection between the smooth behaviour at T = 
and the finite- temperature transition has not been ad- 
dressed before. Our model study also has direct implica- 
tions for the volume-collapse transition of materials such 
as Cen sLan 1 Thn 1 jlfij . and its dependence on magnetic 
field [Fjj ! as explained at the end of this paper. 

We study a generalisation of the periodic Anderson 
Model (PAM) defined by the hamiltonian: 
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In addition to the usual hybridisation and interaction 
terms, it contains a direct hopping between the /- 
orbitals: tff = at. The model reduces to the PAM when 
a = 0. When V = 0, it describes two independent flu- 
ids: free conduction electrons, and a narrow band of /- 
electrons described by the Hubbard model. For simplic- 
ity, our study is restricted to the particle-hole symmetric 
case where both bands are half-filled ((rif) = (n c ) = 1). 
In this case, one has a (renormalised) hybridisation-gap 
insulator when the direct /-/ hopping is small, as studied 
in (for a — and large U this is the Kondo insula- 
tor). For a large enough however, the hybridisation gap 
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closes and the model describes a metal. As shown below, 
the criterion for a metallic ground-state is essentially in- 
dependent of U and reads: a > (V/D c ) 2 (with D c the 
conduction electron bandwith). 

In this article, we study this metallic regime within 
DMFT, focusing on the paramagnetic phase. When V — 
0, the situation is well-documented 0,0 : the f-electrons 
are described by a Hubbard model which undergoes a 
Mott transition. The transition is first-order at finite 
temperature, with a transition line U C (T) ending at a 
critical endpoint (U C ,T C ). The transition line separates 
two different regimes: on one side the /-electrons are 
itinerant while on the other side the /-electrons behave 
as local moments. These two behaviors correspond to two 
locally stable mean-field solutions, which coexist in the 
domain U c i(T) < U < U C 2(T) delimited by two spinodal 
lines. The transition persists down to T = 0: there, the 
quasiparticle weight of the correlated itinerant solution 
vanishes continuously at U C 2(T = 0). The main issue we 
want to address is whether a regime with unscreened local 
moments survives in the presence of a finite hybridisation 
(V 7^ 0) to the broad conduction electron band, and what 
happens to the phase transition. 

DMFT associates to this lattice model a single- 
impurity Anderson model for the /-orbital, subject to 
an effective hybridisation function A c g(iuj n ) which must 
be determined self-consistently. For the case of semi- 
circular densities of states for both the c- and /-electrons 
(corresponding to a large-connectivity Bethe lattice) , this 
self-consistency condition can be written, using the cav- 
ity construction 0, as: 

99,-,/ \ \V — at 2 G c f (iuj n )] 2 
A eS (iu n ) = a 2 t 2 G ff {iu n ) + ^-. ° M n> } 2 

lOJn ~ t 2 G cc (lUJ n ) 

In this expression, G// , G c f and G cc are the different 
components of the on-site interacting Green's function, 
which must be computed self-consistently from the effec- 
tive impurity model. This expression has a transparent 
interpretation: The screening of the /-moment on a given 
site in the local picture of the lattice model has two ori- 
gins reflected in each term of this equation. The first 
term describes the screening due to the motion of the 
/-electrons onto other sites: it is effective only when the 
/-electrons are itinerant, and its vanishing at low-energy 
is associated with the Mott phenomenon. The second- 
term describes the local screening due to the conduction 
electrons. This screening is affected by the /-electron 
motion, resulting in a reduced frequency-dependent ef- 
fective hybridization V e s(ioj n ) — V — od 2 G c f(iuJn)- 

Let us consider first the case T = 0. Physical intu- 
ition suggests that an arbitrarily small hybridisation V 
is enough to screen the local moment through the for- 
mation of a Kondo singlet with the conduction electrons. 
The energy scale associated with screening will be very 
small, but finite, at small V. Hence, at T = 0, the 
hybridization is a singular perturbation when starting 



from the paramagnetic Mott phase, suggesting that the 
T = Mott transition is unstable against the introduc- 
tion of hybridization. This intuition is supported by a 
low- frequency analysis of Eq. (J2J . In a Mott phase with 
unquenched /-moments, the Green's function and self- 
energy behave as: T,f(iuj) ~ I / 'iui , G / / (iui) ~ iui at small 
ui. Inserting this into (|2J), one sees that A e g(iui) ~ iuj 
as u> — > if V = 0, which is consistent with the original 
assumption of a local moment as it implies a gap in the 
hybridisation density of states ImA ff (ui+i0 + ). However, 
as soon as V ^ 0, A e g(iui) tends to a finite (imaginary) 
value as uj — > because of the second term in This 
implies a finite value of ImA e ff (ui + i0 + ) at low- frequency, 
which is inconsistent with a free local moment. Hence, at 
T = and when V ^ 0, the self-energy has a local Fermi- 
liquid form T,f(iuj) ~ iu>(l — 1/Z) + ... for all values of 
U, with Z the /-quasiparticle weight. At large U, Z is 
very small and sets the scale for screening. This yields 
two quasiparticle bands, which read (neglecting lifetime 
effects): 2uj± = (I + aZ)e k ± [(I + aZ) 2 e\ + AZV 2 } 1 / 2 . 
This corresponds to the non-interacting bandstructurc, 
with renormalized parameters: a c g = Za, V c ff — \fZ V. 
It is easily seen that a hybridisation gap is present only 
if V e s > v^efi D c . The quasiparticle weight Z drops out 
from this criterion: hence, the two quasiparticle bands 
overlap and one has a metal when V < yfa. D c , indepen- 
dently of U as announced above. Accordingly, it follows 
from this low-frequency analysis that, at T = 0, the /- 
spectral function is pinned at uj = to its non-interacting 
value: Aff(0) — aT 1 po(V/y/a), for all U, as long as 
< V < y/a D c (with po the non-interacting density 
of states of the conduction band). In this Fermi- liquid 
state, the "large" Fermi surface encounters n c + nt — 2 
electrons per site. This low-frequency analysis can be 
illustrated by a simple calculation using the Gutzwiller 
approximation (GA). In this approach, one optimizes a 
variational energy depending on the probability of double 
occupancy d, and the quasiparticle residue is obtained|l9| 
as Z = 16d(l/2 — d). The results of this approximation 
for our model are displayed in Fig. ^ This figure clearly 
shows that the Brinkman-Rice transition (analogous to 
the U C 2 found in DMFT), at which Z vanishes in the 
Hubbard model (V = 0), is no longer present at finite 
V. We showed analytically that, because V introduces 
a logarithmic singularity (oc Z In Z) in the variational 
energy at small Z , the minimum is always found at a 
finite value of Z. This also allows us to estimate the be- 
haviour of Z at large U U C 2, which has the expected 
exponential suppression characteristic of the Kondo ef- 
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(the prefactor c a y depends 



only weakly on a and V) . We performed a full quantita- 
tive solution of the DMFT equations at T — 0, in order 
to confirm and extend this analysis. For this purpose, we 
used an exact diagonalization (ED) scheme based on the 
Lanczos algorithm and an adaptative discretization of the 
effective bath degrees of freedom [lj . The Green's func- 
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FIG. 1: Comparison of Z vs U for different values of V (from 
the left: V = 0.1,0.2,0.3,0.4) obtained using the GA and 
T — ED (inset). Both methods clearly show that Z never 
vanishes as soon as V ^ 0, in contrast to the V = (Hubbard 
model) case, which displays a transition at the critical value of 
U indicated by the arrows. Note that the GA overestimates Z. 
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FIG. 2: Phase diagram within the IPT approximation. Re- 
gions where localized and itinerant solutions coexist are dis- 
played for (from bottom to top) V=0. 1,0.2,0.3. Only finite- 
temperature transitions exist when 1^/0, resulting in two 
critical endpoints. In contrast, the spinodals of the V = 
case (dashed lines) reach T = at finite critical values of U. 



tion and self-energy obtained from ED (not shown) do 

obey G(iw) iira' 1 p (V/ 'y/a) , S(iw) - - 1/Z) 

at low- frequency for all V =/= 0, from which we obtained 
the quasiparticle weight displayed in Fig. ^ Besides 
confirming the analysis above, we also performed ED 
calculations for increasing and decreasing sweeps in U 
in order to check that no other solution of the DMFT 
equation is present at T = when V ^ 0, besides the 
Fermi-liquid one with screened /-moments and a large 
Fermi surface. This is in contrast to the Hubbard model 
(V = 0) which has a coexistence region \U c i,U C 2\ be- 
tween a Mott-localized and an itinerant solution extend- 
ing down to T = 0. 

While Kondo screening always sets in at T = and 
suppresses the Mott transition, we expect the situation 
to be qualitatively different at finite temperature. In- 
deed, on general grounds, the effect of a perturbation 
(even if singular for the ground-state) is expected to be 
smooth at T ^ 0. As a result, the first-order transition 
and the coexistence region should be robust features of 
the present model as long as V is not too large. Since 
small energy scales are involved (T c is of order ~ D c /40 
for the pure Hubbard model, and the Kondo screening 
scale is tiny at large U), an exact numerical study is 
difficult and we approached the problem using two ap- 
proximate impurity solvers. The first is the iterated per- 
turbation theory (IPT) approximation[2£j , which has 
proven to be semi-quantitatively very successful in the 
study of the Mott transition. IPT is known to overesti- 
mate low-energy scales, and will not be accurate in the 
Kondo regime. The second method is the (dynamical) 
"slave-rotor" (DSR) integral equations [ill], which is able 
to resolve low-energy scales and reproduces the correct 
exponential Kondo scale at large U. The phase diagram 
found within IPT is displayed in Fig. |3 As anticipated, 



the coexistence (hysteretic) region is still present for the 
smaller values of V. AsV increases, its extension is dras- 
tically reduced, and T c decreases, as also found with the 
more accurate DSR solver (Fig.[3]compares the estimates 
of T c in the two methods) . The spectral functions of two 
coexisting solutions are displayed in Fig. one has a 
well-formed Kondo peak corresponding to good screening 
of the local moment while the other one has very small 
(but finite) spectral weight at low energy. In contrast to 
the V = case, we find that the two spinodals no longer 
extend down to T — and that, within IPT, another crit- 
ical endpoint is found at low temperature at which the 
actual first-order transition line terminates. Note that, 
in view of the above analysis at T = 0, the two spinodal 
lines must indeed either end at a lower critical point or 
run away towards infinite coupling. Unfortunately, be- 
cause of the low energy scale involved, we have not been 
able to push the DSR method to low enough tempera- 
tures and firmly establish the existence of a lower critical 
endpoint within this technique. 

We have also studied the effect of a magnetic field on 
the transition between the screened (low-T, itinerant) 
and unscreened (high-T, local moment) regimes. This 
is motivated by recent experiments on the 7 — a transi- 
tion of the Ceo.sLao.iTho.i alloy, showing that the tran- 
sition temperature is decreased by an applied magnetic 
field ■ Dzero et al. [13] pointed out that this can be ra- 
tionalized by approximating the high-T 7- phase as a col- 
lection of almost free localized magnetic moments, while 
assuming that the free-energy of the low-T cv-phase does 
not change appreciably with magnetic field. In Fig. [3)3, 
we display our findings for the temperature associated 
with the local moment spinodal line (lower boundary of 
the coexistence region), as a function of applied field. 
This demonstrates, within the simple microscopic model 
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FIG. 3: Left panel: comparison between spectral densities for 
V (dashed lines) and V = (Hubbard Model, full lines), 
in the screened phase (A) and in the unscreened one (B). 
C: Upper and lower T c of the critical endpoints as functions 
of V within IPT (full lines). Upper T c with DSR (dashed 
line). D. Transition temperature as a function of an applied 
magnetic field H for U=2.7, calculated at T = within ED. 



studied here, that indeed the transition is suppressed by 
a magnetic field as observed experimentally. 

To summarize, we have studied a model which retains 
the key aspects of f-electrons metals, i.e. a narrow cor- 
related band hybridized with a wide uncorrelated con- 
duction band. We found a first-order transition at fi- 
nite temperature between a screened phase and a local- 
moment phase, ending in two critical endpoints. We 
demonstrated that, at zero-temperature, the hybridiza- 
tion is a relevant perturbation, so that the Mott transi- 
tion is suppressed by Kondo screening. This is qualita- 
tively consistent with the KVC picture. These findings 
can be put in the broader context of the orbital-selective 
Mott transition (OSMT) , which attracted a lot of at- 
tention recently [23ll24|. In a general two-band case, the 
self-energy 2 takes a matrix form and the Mott transition 
is signaled, when approached from the metallic side, by a 
low-frequency singularity in to. I— = Z~ 1 tu+- ■ ■ . An 
OSMT is characterized by Z having one zero-eigenvalue 
while the other one remains finite. We emphasize that 
this is a basis-independent notion. In our simple model, 
Z is diagonal, with Z cc — 1, and hence an OSMT phase 
corresponds simply to the vanishing of Zft. At T = 0, 
this does occur for zero-hybridisation, but Kondo screen- 
ing prevents this to happen when a finite hybridization 
is turned on. 

Finally, we comment on the qualitative relevance of 
our results for f-electron materials. There, the contribu- 
tion (F e ) to the free-energy from the electronic degrees 
of freedom that are active close to the transition have 
to be added to the contributions from all other bands 
and ions, which can be approximated by an elastic con- 
tribution (see e.g 0). As a result, the volume-collapse 
transition does not correspond to a true divergence of 



the response function \ e — —d 2 F e /dv 2 of the active elec- 
tronic degrees of freedom (with v the unit-cell volume). 
Rather, it will take place[25j when \e = B q /vq, with 
Bq and v a typical bulk modulus and unit-cell volume. 
Hence the critical temperatures of the upper and lower 
endpoints are shifted upwards and downwards, respec- 
tively, by elastic terms. We therefore conclude that there 
are two generic situations that are consistent with our 
results: either the material displays a first-order transi- 
tion down to T = (for softer materials, with smaller 
values of Bq) or it will display two critical-endpoints (for 
harder materials) (see also 0). Experimental studies (ltij 
suggest that the latter case may be realized in the Cerium 
alloy Ceo.8Lao.iTh .i, in which alloying acts as a "nega- 
tive pressure" , thus allowing for an investigation of the 
localization-delocalization transition at lower tempera- 
tures than in pure Cerium. We emphasize that our re- 
sults also imply that a T = quantum-critical valence 
transition is a non-generic case that requires the tuning 
of an extra parameter. 
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